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' The concept of environmental interface is defined and analyzed from the point of view of the 

Pj^ , possible source of non-standard behaviour. The energy balance equation is written for the interface 

where all kinds of energy transfer occur. It is shown that under certain conditions, the discrete 
version of the equation for the temperature time rate turns in to the well-known logistic equation and 
the conditions for chaotic behaviour are studied. They are determined by the Lyapunov exponent. 
The realistic situation when the coefficients of the equation vary with time, is studied for the Earth- 
environment general system. 
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I. INTRODUCTION 



The concept of "balance", either global or local in any given context, it is undoubtedly the cornerstone in the 
increasing number of environmental problems. The question is: Why are the environmental problems in the focus 
' now? One particular answer can be found in a hierarchy of the main scientific problems in this century. According to 
I/"") , the physicists, the world scientific community, will be occupied in this century among others, by the environmental 
problems that are primarily expressed through the problem of climate changes [1, 2] as well as many other linked 
£^ problems that comprises a wide range of time and spatial scales. This is the first time in the history of science that 
the environmental problems take the place at the research front of the sciences. The question, why it is happening 
now and why it will go on happening in the future, could be answered by the well known fact that in the scientific 
. as well as in other worlds the main "dramatic event" takes place at the interface between either two media or 
two states [3]. The field of environmental sciences is abundant with various interfaces and it is the right place for 
application of new fundamental approaches leading towards better understanding of environmental phenomena. We 
define the environmental interface as an interface between the two, either abiotic or biotic environments each, which 
are in a relative motion exchanging energy through biophysical and chemical processes and fluctuating temporally 
and spatially regardless of its space and time scale, fn our opinion, this definition broadly covers the unavoidable 
multidisciplinary approach in environmental sciences and also includes the traditional approaches in sciences that 
are dealing with the environmental space less complex than any one met in reality. The wealth and complexity of 
processes at this interface determine that the scientists, as it often seems, are more interested in a possibility of 
non-linear " dislocations" and surprises in the behaviour of the environment than in a smooth extrapolation of current 
trends and a use of the approaches close to the linear physics [3-6]. To overcome the current situation we have 
to do the following: (a) establish a way in approaching the non-linear physics and the non-linearity in describing 
the phenomena in physics as well as environmental sciences, and (b) solve or, at least, understand the problem of 
predictability. These two problems are problems par excellence of the methodology. Their successful solving will help 
us avoid the current problems in mathematical, physical, biological and chemical interpretation of the nature, fn the 
next section, we write down the equation for the energy balance at the environmental interface. We show that under 
certain conditions, its discrete version becomes well-known logistic equation leading to the chaotic behaviour of the 
temperature. The conditions for such behaviour arc investigated and particular examples studied. The analysis is 
based on the Lyapunov exponent. 
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II. ENERGY BALANCE EQUATION FOR ENVIRONMENTAL INTERFACE 

We see the outside world, i.e. the world of phenomena (ambience) from the observer's perspective (its inner world). 
In the ambience there are systems of different levels of complexity and their environments. System in the ambience 
is a collection of precepts while whatever lies outside, like the component of a set, constitutes the environment [7,8]. 
The "fate" of science lies in the fact that it is focused on the system [8]. Furthermore, to be able to anticipate 
something, we describe the system by the states (determined by observations) while the environment is characterized 
through its effects on the system. The environmental interface as a complex system is a suitable area for occurrence 
of the irregularities in temporal variation of some physical or biological quantities describing their interaction. For 
example, such interface can be placed between: human or animal bodies and surrounding air, aquatic species and 
water and air around them, natural or artificially built surfaces (vegetation, ice, snow, barren soil, water, urban 
communities) and atmosphere [4-6,9,10], etc. The environmental interface of different media was recently considered 
for different purposes [3-6, 11-13]. In these systems visible radiation provides almost all of the energy received on the 
environmental interface. Some of the radiant energy is reflected back to space. The interface also radiates, in the 
thermal waveband, some of the energy received from the sun. The quantity of the radiant energy remaining on the 
environmental interface is the net radiation (the net radiation energy available on the surface when all inward and 
outward streams of radiation have been considered as seen in Fig. 1 where all kinds of energy are expressed in terms 
of the flux density), which drives certain physical processes important to us. The energy balance equation may be 
written as 

c~^=R-H-E-S (2.1) 

where Ci is the environmental interface soil heat capacity per unit area, Ti is the environmental interface temperature 
(EI), t is the time,i? and E are the sensible and the latent heat, respectively, transferred by convection, and S the 
heat transferred by conduction into deeper layers of underlying matter. 

Figure 1 

Schematic diagram of terms included in the energy balance equation for environmental interface: (a) net radiation 
(R) that includes (1) visible radiation, (2) infrared radiation of underlying matter and (3) infrared counter radiation 
of the gas; (b) sensible(ff) and latent heat (E) transferred by convection and (c) heat (S) transferred by conduction. 

The sensible heat is calculated as Ch^Ti — T a ) where Ch is the sensible heat transfer coefficient and T a (t) is the gas 
temperature given as the upper boundary condition. The heat transferred into underlying soil material is calculated 
as CoiTi—Td) where Cd is the heat conduction coefficient while Td(t) is the temperature of deeper layer of underlying 
matter that is given as the lower boundary condition. Following Bhumralkar [14] the net radiation term in Eq. (1) 
can be represented as Cn(Ti — T a ) where Cr is the radiation coefficient. According to [15], for small differences of T a 
and Td, the expression for the latent heat can be written in the form Cl f{T a ) [b(T[ —T a ) + b 2 (Ti — T a ) 2 }/2. Here C^is 
the latent heat transfer coefficient, f(T a ) is the gas vapor pressure at saturation and b is a constant characteristic for 
a particular gas. Calculation of time dependent coefficients , and can be found in [16]. After collecting all terms in 
Eq. (1) we get 

^ = Ai(Ti - T a ) - A 2 (T l - T a f - A 3 (T t - T d ) (2.2) 
at 

where A\ = [Cr — Ch — CLbf(T a )]/ci, A 2 = b 2 f (T a ) / (2ci) , and A3 = Co/ci- This is a non-linear Riccati type 
differential equation that practically always has to be solved numerically, i.e., 

DT, = F„ (2.3) 

(3) where D is the finite difference operator defined as DT, = (T^ n+ i — Ti >n )/T>t, n the time level, Tit is the time step 
and F„ is the r.h.s. of Eq. (2.2) defined at the nth time level. 

III. NUMERICAL RESULTS ON CHAOS IN TEMPERATURE ON ENVIRONMENTAL INTERFACE 



We consider Eq. (2.2) with the lower boundary condition, at some time interval, given in the form 

T d = T a -(c,/C D )^ 



(3.1) 
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expressing slow temperature changes in both the environment and underlying material. Then equation becomes 



§ = A ti - A 2 e (3.2) 

where £ = Tj — Ta, and Aq = A\ — C'o/ci. Substituting the time derivative by the finite difference operator D in this 
equation and after some transformations, we obtain 



r n+1 = A?r„(i-r n ) (3.3) 

where V = (A^/A\)^, A\ = 1 + A T)t and A v 2 — A 2 T>t. This equation has the same form as the well known logistic 
difference equation r n+ i = f3T n (l — T n ) where (3 is a constant (see, e.g. [17] , among others). Conditions for occurrence 
of deterministic chaos for logistic mapping given by this equation are: (a)0 < V < 1 and (b) 3.57 < (3 < 4. Here we 
have to bear in mind that A\ depends on discrete "time" n. However, chaos is still expected, although we can give 
no quantitative prediction of the parameter regime. We analyze Eq. (3.3) in the following way. With TDt p = I/Aq we 
indicate the scaling time range of energy exchange at the environmental interface including coefficients, that express 
all kind of energy reaching and departing the environmental interface. For any chosen time interval, for solving Eq. 
(3.3), there always exists T)t p j = Min\Dt p {c.i, Cr, Ch, Cl)] when energy at the environmental interface is exchanged 
in the fastest way by radiation, convection and conduction. If we define dimcnsionless time r = T)t/T)t p j, then Eq. 
(3.3) becomes 

r„ +1 = (i + r)r„(i-r„). (3.4) 

Regarding to order of magnitude of parameters included in coefficients A\ and A%, the condition (a) is always 
satisfied. From (b) we get the interval ((3 — 1 + r ) where the solution is chaotic, i.e., 2.57 < r < 3. We analyze now 
the occurrence of the chaos in solution of Eq. (3.4). Since conditions (a) and (b) are necessary but not sufficient for 
identification of the chaos, we shall calculate values of the Lyapunov exponent for situations when these two conditions 
are satisfied. Here we define the Lyapunov exponent Xl for the case of Eq. (3.4), which has a single degree of freedom 
r which depends on discrete "time" j following the form [18] 

1 - 

X L = Um n ^ 00 -yin\P'{r j )\ (3.5) 
i=i 

where P(T) = (1 + r) T(l - L). 

To examine how changes in Aq and A 2 determine the irregularities in behaviour of the temperature at the 
environmental interface obtained from the energy balance equation, we perform numerical experiment. We solve Eq. 
(3.2) for different values of Aq and A 2 , varying them in a broad range of energy exchange coefficients. After that 
we calculate the corresponding Lyapunov exponent Xl as a function of r. Figure 2 depicts dependence of Lyapunov 
exponent Xl on t. It takes mostly positive values indicating chaotic fluctuations of T. However, inside of the chaotic 
interval there are a lot opened periodical "windows" where Xl < 0. 

Figure 2 

Dependence of Lyapunov exponent Xl on the dimcnsionless time r. 

Now, we give a specific example of the above analysis applied to the interface of the Earth and its environment. It 
is now accepted that the Earth is a complex system which consists of the biota and their environment. These two 
elements of the system are closely coupled: the biota regulate the environment (e.g., climate on planetary scale) 
and, in turn, the environment restricts the evolution of the biota and dictates what type of life can exist. Changes 
in one part will influence the other, being opposed by negative feedback or enhanced by positive feedback, and this 
may lead to oscillation or chaos in the system. In the system considered, the values of the dimcnsionless time r 
for long term atmospheric integration are in the interval 1 < t < 3. Corresponding bifurcation map given in Fig. 
3 depicts the regions with the chaotic fluctuations of temperature at the environmental interface biota-surrounding air. 

Figure 3 

Bifurcation diagram of T as a function of dimcnsionless time r in the long term atmospheric integration. 
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IV. CONCLUSION 



The main aim of the report was to indicate the possibility of chaotic behaviour of the temperature at the interface 
of two media where the energy is exchanged by all three known mechanisms. The conditions for such phenomenon 
are discussed. Once the dimensionless time is introduced, one can study the problems on rather different scales in 
the same manner (scaling approach). 
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